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Алгоритм решения задачи 
об оптимальньх назначениях методом границ и ветвей 


Наведено нову (комбінаторну) модель задачі про призначення. Досліджено особливості використання 
методу гілок та меж для розв'язування задачі про призначення. Поліпшено оцінку допустимих множин у 
методі гілок та меж; розроблено алгоритм розв'язування задачі та проілюстровано його на прикладі. 
Ключові слова: метод гілок та меж, задача про призначення, оцінювання підмножин, 
оптимізація на перестановках. 


А пем (сопабіпаюогіаї) подє! аззідптепі ргобіет. ТПе Кеагигез об "ре теїрод Бгапсі апа бошпа бог 50Їміпе 
Фе аззієптепі ргобіет. Ппаргоуей аз5е55тепі ої айтіз5ібіє 5еї5 іп Бгапсб апа Боппа, Ше аїсогіїйт Їог 
50ЇІуїпє, Фе ргобіет апа Шибвігаїе її м'їїб ап ехатріе. 

Кеу ууогая: теїпод Бгапсі апа боицпа, Фе (азКк ої арроіпітепі, еуаїцайоп ої 5иБееїз, оріітігацйоп 
оп регтиїайоп5. 


Показана полная (комбинаторная) модель задачи о назначениях. Исследованьт особенности использования 
метода ветвей и границ для решения задачи о назначениях. Улучшена оценка допустимьїх множеств в 
методе ветвей и границ; разработан и проиллюстрирован алгоритм решения задачи на примере. 
Ключевьге слова: метод ветвей и границ, задача о назначениях, оценка подмножеств, 
оптимизация на перестановках. 


Вступ 
Задача про призначення - одна з класичних задач цілочислової оптимізації | 1). 
Її суть полягає в тому, що є и видів робіт та п кандидатів для їх виконання (виконавців). 
Вважається, що кожен з кандидатів іє ..,, - (1,....п) може виконувати будь-яку 
роботу іє, з В. , при цьому с, - ефективність виконаної роботи /-го виду і-м 


кандидатом. 

Необхідно так розподілити кандидатів на виконання робіт, щоб кожен з канди- 
датів одержав єдине призначення, кожна з робіт одержала єдиного виконавця 1 су- 
марна ефективність, пов'язана з призначеннями, була максимальною. 
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Мета даної статті - розробити нову (комбінаторну) модель задачі про призна- 
чення. 

Розробка алгоритму за методом гілок та меж є досить актуальною для задач 
такого типу, оскільки дає можливість використання його у розв'язуванні більш склад- 
них задач (наприклад, задача теорії розкладів) |21. 

Побудуємо нову модель задачі в рамках евклідової комбінаторної оптимізації (3-5. 

Введемо змінні: 


1, якщо і- та робота виконана ) -- м виконавцем; 


х 
7. 0, в іншому випадку. 


Цільова функція має вигляд 


к9-У У гу, -» тах. (1) 


із у-1 
Повинні виконуватись такі умови: 
- на кожну роботу призначається тільки на один виконавець 


Уску нм ен (2) 
і«і1 

- кожен виконавець призначається тільки на одну роботу 
У ан МРе І (3) 


1-1 
я. і 1011) й 
Зауважимо, що остання умова може бути записана комбінаторно. 


Розглянемо мультимножину С ог) з основою 5 (с) - (0,1) та первин- 


ною специфікацією |С|- (п -пу) 
Тоді допустимий розв'язок задачі про призначення є впорядкованою /-вибіркою з 
множини С, де А - п", тобто елементом загальної множини переставлень Е, , (с) : 
хе( худ бнь Ки ЛЕ Воз (С): (4) 


Задача (1) - (4) є моделлю задачі про призначення як лінійної умовної повністю 
комбінаторної задачі на перестановках |21. 


Оцінювання допустимих підмножин 
в методі гілок та меж 


Моделі задач типу задачі про призначення, але з додатковими умовами (задачі 
теорії розкладів |2|, є актуальними, але до них не можна застосовувати поліноміальні 
алгоритми (угорський метод |6|) тощо). Тому вбачається для таких задач (почавши з 
задачі про призначення) дослідити алгоритм у класі методу гілок та меж |4-61. 

Метод гілок та меж (МГМ) - універсальний метод, який можна застосувати до 
будь-якої оптимізаційної задачі. У випадку задачі про призначення з матрицею 


ефективності С - (с,) будемо використовувати модель (1) - (4) цієї задачі у вигляді 


лінійної умовної задачі комбінаторної оптимізації на переставленнях. 
У подальшому будемо вважати, що елементи С впорядковані по незростанню: 


ор Р Р (5) 
тобто 8) 2-2, 2-78, 75 Є Зо ТТ 810. 
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Суттєвим для ефективності МГМ є оцінювання допустимих підмножин. 

Як відомо |4-6|, в задачі максимізації на множині ) функції К (х) , оцінка для 
підмножини Р, с ) - це число у, 2 К(х) Ухе ), У задачі мінімізації знак не- 
рівності міняється на протилежний. 

Для організації оптимізації за методом гілок та меж необхідно визначити: 
1) спосіб оцінювання допоміжної підмножини Ї) ; 2) спосіб утворення 0); 3) прави- 
ло (правила) відсікання безперспективних чи порожніх підмножин Щ). 

В рамках цієї роботи розглянемо оцінювання, а при виконанні алгоритму утво- 
рення та відсікання 0, як це розглядалося в |4-6|. 

Оцінювання підмножини 0), можемо здійснювати за допомогою підходу ви- 
кладеного в теоремі 5.3 |4, с. 129| у вигляді 

дауче, 

де у, - це сума частини доданків у цільовій функції Б (х) (1), що визначена 
вже заданими невідомими, які визначають Д,; с - оцінка невизначеної частини 
цільової функції КГ (х) для Ф,, тобто тих її доданків, що містять незадані в Д, змінні. 

Оцінка 2, підмножини Д, може бути покращена за рахунок покращення с! 
внаслідок викреслення з матриці С елементів, які стоять в рядках і,...,і, та стовпцях 
ЕРИ ДУ 

Матрицю, утворену з матриці С - С (і, - /, 2 0), позначимо СХ. Вона утво- 
рюється викресленням з С"! рядка і,, стовпця /,. Підмножина множини допустимих 


розв'язків (2)- (4) 1" задається значеннями змінних х, зх, з..тх,, ЗІ 


йЛ ьь СТОЛ 
Елементи матриці чі перепозначимо і перенумеруємо так: 
сеовавО (6) 
ПОРО а» 
не вказуючи залежності від пар (і 1); т - 1,2,....Ї, якщо це не викликає плута- 


нини. 


Для зручності елементи матриці С" будемо позначати як упорядковану множи- 


І 
п-І,п- І 


І . . 
ну С ар . При необхідності будемо використовувати елементи з С" з 


індексами з С, тобто с,, з Су: 
о | з . | . . . . -. . 
Нехай І Й 94 -і нан у Покращення оцінки допустимої під- 
множини сформульовано у вигляді такого твердження. 
Теорема 1 (про оцінку в задачі максимізації). Оцінкою 2/7," підмножини 


Р723/ може слугувати величина 


ці аЙ 


ло ума є, ; (7) 


й й 


І 
зе | ЛллнЛ ЛллеЛ у 
де С Пи ря а Удай обчислюється так Урі ш Сг, . 


1ч1 їч1 


16 «Искусственньтй интеллект» 2013 Хо2 


Алгоритм розв зування задачі про оптимальні призначення методом гілок та меж | л 


Доведення. Якщо х,, 2 Х,, 2-7 Х,, - 1 визначають підмножину 277", то 


2/1 


для МУ є ро «Л 


по 


удо Л 
АХ х)з Уа АЛ кох, і 


іє! 
є! 


при виконанні умов (2) - (4). Що, очевидно, дає нерівність Ра суХу 2 С,, ТОБТО 


іє!! 
у! 


ло чо жд Лін : 
Е(хуєирйлт зар аби Мхе ДИ, Що і треба було довести. 


цій ЛИЙ ЛІ ці 
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Крок 1. Упорядковуємо та пронумеруємо одним індексом по незростанню еле- 
менти матриці ефективності С: 


Що 
ІМ 
М 


(8) 


де с, є они 

Крок 2. За допомогою «жадібного» алгоритму знаходимо допустимий розв'я- 
зок задачі про призначення і відповідне значення сумарної ефективності як суму з п 
елементів матриці ефективності ЕК. Цей алгоритм може бути таким. Нехай на по- 


чатку 5, -0. 

Знаходимо найбільший елемент матриці, додаємо його до Б, викреслюємо рядок 
і стовпець, в якому знаходиться знайдений елемент. З отриманою матрицею робимо 
те, що й з вихідною, поки не залишиться один елемент, який також додаємо до ЕК. 
РЕ, відіграє роль поточного рекордного значення цільової функції задачі, тобто 


найбільше зі знайдених значень. 
Оскільки знайдене ЕК, не є гарантовано оптимальним, то його можна спро- 


бувати покращити у наступний спосіб. Спираючись на той факт, що в оптимальному 
розв'язку містяться доданки с, , які не менші, ніж и найбільших елементів в упоряд- 
куванні (8), починати визначення Б, можна не тільки з найбільшого, а з будь-якого з 
п «найбільших» елементів в (8), вибираючи в якості Б, з одержаних и сум найбільшу. 

Крок 3. Задаємо початкові значення /-0, т -0, / - номер рівня дерева, т - 
номер елемента матриці, що використовується в підмножині. Позначаємо для одно- 
манітності формул С - С", р «Фо. 

Зауваження. В кожній підмножині допустимих розв'язків, що розглядаються в 
алгоритмі, значення г своє, тому для 1)? т позначимо т, (далі і,, /, при наявності 
і, д опускаємо). 

Крок 4. Збільшуємо / на одиницю / :-/-1Ї. 

Крок 5. Збільшуємо на одиницю номер елемента множини С" п Ре ЗР 


Якщо т - (п-І) мІє, та І - п, то перехід на крок 12; при / «и хоча б для 


- 5 2 
одного / - перехід на крок 11, інакше (тобто коли т//// є (п-Г) ) -крок б. 
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Крок 6. Коефіцієнту цільової функції ставимо у відповідність елемент множи- 
ни С (одиницю) в якості відповідного значення х,. Для цього вибираємо наступний 
елемент с"! матриці С"! в порядку (6), який в С має позначення с, - Задаємо х,, 21. 
Виокремлюємо підмножину рі б з підмножини ри за Далі таку підмножину 
позначатимемо Щ0/7". Утворюємо матрицю С" з С"!, викресливши рядок і, 
стовпець /,, (тобто переходимо до наступного «С -- 1) -го) рівня дерева галуження). 


Елементи С" перепозначимо згідно з (б). 
Крок 7. Якщо / я п, то перехід на крок 8. Інакше (п - ) , обчислюємо Б (х) 


шо - з і ке » За і 
при ху ЗХ, ЗХ, 71» інші х, 2-0. Якщо Е(х)2 Ю, то Б, - К(х). Перехід на 
крок 11. 

Крок 8. Знаходимо для підмножини 19/77" оцінку 2, що обчислюється за 


формулою (7). 


Крок 9. Порівнюємо оцінку 2.77" 


-й 


і ЛІ 
із значенням ЕК. Якщо 2; 2ЙЕ,, пере- 


ходимо на крок 4, інакше на крок 19. 
Крок 10. Оскільки ві Уг «Є Б, то подальше галуження поточної підмножини 


ДЛ припиняється, вона відсікається. 


ій 

Крок 11. Якщо г/// - (п -1І У» , то перехід на крок 12, інакше на крок 5. 

Крок 12. Повернення на (1-1) -Й рівень дерева. Для цього використовуємо 
останню розглянуту підмножину цього рівня 0/7, "!. Тобто зменшуємо / на одини- 
цю: / :- 1 -1. Перехід на крок 5. 

Крок 13. Вивід відповіді - значення цільової функції Б, та набір векторів х, 


множина всіх оптимальних переставлень. Зупинка. 

Проілюструємо елементи виконання алгоритму на прикладі. Нехай задана 
матриця ефективності С для задачі про призначення. Знайти розв'язок задачі за 
допомогою МГМ. 


Са 


что Кк в ча 
- со с 10 
ь о ь с 


3 
3 
2 
4 


Упорядковуємо по незростанню елементи матриці ефективності С: 
9-сь»ТаесьсТасу»ббаСсзвВаСся»У ВО ВУва»фНо З4ас,2 
з4-с, зас, »ЗасуаЗасуаЗасц»2с 22 со 
За допомогою жадібного алгоритму знаходимо 5, починаючи визначення з п 
найбільших значень су» В, -21. 
ПОРО УРаК ОК СЛЬКУТ (9) 
На першому рівні галуження привласнює х;; - 1, утворюємо ДР. Знаходимо для 


р? перший доданок оцінки (7) у; - 9 та оцінку 2 - 22.» Е,. Продовжуємо галуження, 
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переходимо до другого рівня (/ - 2), утворюємо Ді: ху, -1, у; з14, 2, -212К. 


. бо б М 324. і 324 || 334 
Аналогічно на третьому рівні (/ - 3) галуження - Ді, 1 Х4 - 1, Уцо - 17, 2» -20«К,, 


тому підмножина 7 відсікається (на рис. І знак 9 ). Переходимо до наступного 


елемента матриці на тому ж рівні (1 б 3) і так далі (рис. 1). 


Ду «98 е32»к ру зі ат» к, ру 18,28 -20«Е, 


142 14 4 142 2 /142 142 


о (0) Хз -1 Ретро! » 


Рисунок 1 - Фрагмент дерева розв'язків для трьох підмножин 


Для наведеного прикладу МГМ показав, що покращення Й, не відбулося. Отже 


2 Ре - аг ен рен, 
розв'язком є БЕ, - 21 при призначенні, що задається формулою (9), тобто хі, в Хо З Ху З Ху: 


Висновки 


У статті розроблено один з алгоритмів пошуку розв'язків методом гілок та меж 
для задачі про призначення. Побудовано нову математичну модель задачі. Сформульо- 
вано та доведено твердження про покращення оцінки допоміжних підмножин. Детально 
описано способи галуження та відсікання порожніх та безперспективних підмножин. 

В подальшому доцільно розглянути поширення способів галуження, оцінюван- 
ня і відсікання, наведених в роботі, до задач теорії розкладів |21. 
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ВКЕ5СМЕ 
М.Х. Іеопоуа 
Асогійт Зоіуїпя їле Ргобіет о/ Оріїта! Кипсіоп Бу Метоа 
ОЇ ВкапсП апа Воцпа 


ТРе агіїсіе деаіз муїії Фе арріїсанйоп ої опе ої Ше по5і уег5аціе птеШод5 Гог 50Їуїпе 
орійтігайоп ргобіетя -- пеїод бгапсб апа бопла. Модеї ргобіетз5 ої їБі8 гуре, Би мір 
аддійопа! сопаїйоп5 (га5к 5спедийпо їБеогу) аге геіеуапі, Бик Шеу сап пої Бе цей 
роїупотіа! аїдогіртя8. Тре 5киду ої 5исВ ргобіете5 аге геЇсуапі аЇсогіїрт5 їп Фе сіає5 об 
ргапсрез апа Боцпаз. а соплесіїоп угіїб їБі5 аЇдогійрт бгапсб апа боцла тефодй їог Фе 


азвідптепі ргобіегп. паргоуедй езійтабоп ої зесопдагу зибзеї ДО, , дебпед апхіагу зибееї 


еуашайоп теїрод, Фе ухау ої ії5 езгабіяптепі апа гиіе ргипіпє ппрготізіпо ог егаріу 
зибясів Д),. 


Стаття надійшла до редакції 09.04.2013. 
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